Scaling properties of azimuthal anisotropy 
in Au+Au and Cu+Cu collisions at ^/snn = 200 GeV 
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Detailed differential measurements of the elliptic flow for particles produced in Au+Au and 
Cu+Cu collisions at ^/snn ~ 200 GeV are presented. Predictions from perfect fluid hydrodynamics 
for the scaling of the elliptic flow coefficient V2 with eccentricity, system size and transverse energy 
are tested and validated. For transverse kinetic energies KEt = mr — m up to ~ 1 GeV, scaling 
compatible with the hydrodynamic expansion of a thermalized fluid is observed for all produced 
particles. For large values of KEt, the mesons and baryons scale separately. A universal scaling 
for the flow of both mesons and baryons is observed for the full transverse kinetic energy range of 
the data when quark number scaling is employed. In both cases the scaling is more pronounced in 
terms of KEt rather than transverse momentum. 

PACS numbers: PACS numbers: 25.75.Dw 



Quantum Chromodynamics calculations performed on 
the lattice (LQCD) indicate a transition from a low- 
temperature phase of nuclear matter, dominated by 
hadrons, into a high-temperature plasma phase of quarks 
and gluons (QGP) (Ij. For matter with zero net baryon 
density, this phase transition has been predicted to oc- 
cur at an energy density of ~ 1 GeV/fm"^ or for a critical 
temperature Tc ~ 170 MeV Recent estimates from 
transverse energy {Et) measurements at the relativistic 
heavy ion collider (RHIC) have indicated energy densities 
of at least 5.4 GeV/fm^ in central Au-|-Au collisions 
Thus, an important prerequisite for QGP production is 
readily fulfilled at RHIC. Indeed, there is much evidence 
that thermalized nuclear matter has been created at un- 
precedented energy densities in heavy ion collisions at 
RHIC iiiilliiEl 

Hydrodynamics provides a link between the funda- 
mental properties of this matter (its equation of state 
(EOS) and transport coefficients) and the flow pat- 
terns evidenced in the measured hadron spectra and az- 
imuthal anisotropy 0,0,0,0,0]. Experimentally, 
such a momentum anisotropy is commonly characterized 
at mid-rapidity, by the even order Fourier coefficients 

0113, 

I'n^ (e^^*--*"-)), n = 2,4,.., (1) 

where 0^ represents the azimuthal emission angle of a 
particle, ^jip is the azimuth of the reaction plane and 
the brackets denote statistical averaging over particles 
and events. 

At low transverse momentum (pr ^ 2.0 GeV/c) the 
magnitude and trends of elliptic flow, measured by the 
second Fourier coefficient V2, is found to be under- 
predicted by a hadronic cascade model 0|. By contrast, 
a broad selection of the data showed good quantitative 
agreement with perfect fluid (very low ratio of viscosity 



to entropy) hydrodynamics 0,0,0,0] and a transport 
model calculation which incorporates extremely large 
opacities jlij . For higher pT, quark coalescence from 
a thermalized state of flowing partonic matter l20j has 
been found to be consistent with the data 0, l22l ]. To- 
gether, these results provide evidence for the production 
of a strongly interacting QGP whose subsequent evolu- 
tion is similar to that of a "perfect" fluid 0>S0,0]- 

Systematic theoretical and experimental studies of the 
influence of model parameters are now required to gain 
more quantitative insight on the transport coefficients 
and the EOS for this strongly interacting matter. The 
range of validity of perfect fluid hydrodynamics is af- 
fected by the degree of thermalization 23] and the onset 
of dissipative effects ^^,.24, 25]. These questions can be 
addressed by investigating several scaling predictions of 
perfect fluid hydrodynamics 001210 ■ 

In the hydrodynamic model, elliptic flow can result 
from pressure gradients due to the initial spatial asym- 
metry or eccentricity e = ((y^ — a;^))/((y^ -I- a;^)), of the 
high energy density matter in the collision zone. The ini- 
tial entropy density S{x,y), can be used to perform an 
average over the x and y coordinates of the matter in the 
plane perpendicular to the collision axis. Here, x points 
along the impact vector and y is orthogonal to x. For a 
system of transverse size R {1/R — \Jl/{x'^) + \/{y'^)), 
this flow develops over a time scale ^ R/ {cg) for mat- 
ter with an average speed of sound Cg. Thus, the initial 
energy density controls how much flow develops globally, 
while the detailed development of the flow patterns are 
largely controlled by e and Cs . 

An important prediction of perfect fluid hydrodynam- 
ics is that the relatively "complicated" dependence of 
azimuthal anisotropy on centrality, transverse momen- 
tum, rapidity, particle type, higher harmonics, etc can 
be scaled to a single function [iH 0] . Immediate con- 



4 



sequences of this 0, are that: (i) V2 scal- 

ing should hold for a broad range of impact parameters 
for which the eccentricity varies, i.e. V2{pt)/^ should 
be independent of centrality; (ii) V2(j>t) should be inde- 
pendent of colliding system size for a given eccentricity; 
and (iii) for different particle species, V2{KEt) at mid- 
rapidity should scale with the transverse kinetic energy 
KEt = TTiT — TTL, where rriT is the transverse mass of the 
particle. 

We use high statistics V2 data to test these scaling pre- 
dictions and explore constraints for the range of validity 
of perfect fluid hydrodynamics. The measurements were 
made at ^/sWn = 200 GeV with the PHENIX detec- 
tor at RHIC. Approximately 6.5 x 10** Au-|-Au and 
and 8.0 x lO'' Cu+Cu minimum-bias collisions were an- 
alyzed from the 2004 and 2005 running periods, respec- 
tively. The collision vertex z, along the beam direction 
was constrained to be within \z\ < 30 cm. The event 
centrality for Au-|-Au collisions was determined via cuts 
in the space of Beam-Beam Counter (BBC) versus Zero 
Degree Calorimeter analog response jsj]. For Cu+Cu 
only the amplitude of the BBC analog response was used. 
Charged hadrons were detected in the two central arms 
{\v\ ^ 0.35). Track reconstruction was accomplished us- 
ing the drift chambers and two layers of multi-wire pro- 
portional chambers with pad readout (PCI and PC3) 
located at radii of 2 m, 2.5 m and 5 m, respectively [30j . 

The time-of-flight (TOF) detector positioned at a ra- 
dial distance of 5.06 m, was used to identify pions (tt*), 
kaons {K"^) and (anti)protons {p)p. The BBCs and TOF- 
scintillators provided the global start and stop signals. 
These measurements were used in conjunction with the 
measured momentum and flight-path length to generate a 
mass-squared distribution |32l | . A momentum dependent 
±2a cut about each peak in this distribution was used 
to identify tt^, and {p)p in the range 0.2 < px < 2.5 
GeV/c, 0.2 < pT < 2.5 GeV/c and 0.5 < pr < 4.5 
GeV/c, respectively. A track confirmation hit within a 
2.5(7 matching window in PC3/T0F served to eliminate 
most albedo, conversions, and resonance decays. 

The differential elliptic flow measurements for charged 
hadrons and identified particles were obtained with the 
reaction plane method. This technique correlates the az- 
imuthal angles of charged tracks with the azimuth of the 
event plane $2, determined via hits in the two BBCs po- 
sitioned symmetrically along the beam line, covering the 
pseudo-rapidity range 3 < \t]\ < 3.9 jl^. A large 77 gap 
between the central arms and the particles used for reac- 
tion plane determination reduces the influence of possible 
non-flow contributions, especially those from jets. Values 
of V2 were calculated via the expression 



V2 



(cos(2(</.p - ^2))) 
(cos(2($2 - $i?p))) 




(2) 



where the denominator represents a resolution factor that 
corrects for the difference between the estimated $2 and 



Pj (GeV/c) 

FIG. 1: V2 vs. pt for charged hadrons obtained in (a) Au+Au 
and (b) Cu+Cu collisions for the centralities indicated, (c) 
V2{centrality,pT) divided by k—3.1 (see text) times the pr- 
integrated value V2(centrality) for Au+Au and Cu+Cu. 



the true azimuth ^b.p of the reaction plane |2ll |33j. 
The estimated resolution factor of the combined reaction 
plane from both BBCs 21] has an average of 0.33 (0.16) 
over centrality with a maximmn of about 0.42 (0.19) for 
Au+Au (Cu+Cu). The estimated correction factor for 
the V2 measurements (i.e. the inverse of the resolution 
factor) ranges from 2.4 (5.5) to 5.0 (13). Relative system- 
atic errors for these V2 values are estimated to be ~ 5% 
and ~ 10% for Au+Au and Cu+Cu, respectively. 

Figure ^ shows the differential V2{pt) for charged 
hadrons obtained in Au+Au and Cu+Cu collisions. The 
V2{pt) results exhibit the familiar increase as collisions 
become more peripheral and the pr increase |^ ^ 1^. 
We test these data for eccentricity scaling by dividing 
the differential values shown in Fig. ^ by the V2 inte- 
grated over the px range 0.3-2.5 GeV/c for each of the 
indicated centrality selections. The hydrodynamic model 
predicts that this ratio is constant with centrality and in- 
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FIG. 2: (a) «2 vs pr and (b) 112 vs KEt for identified particle 
species obtained in minimum bias Au+Au collisions. The 
STAR data are from Refs. 



FIG. 3: (a) v^juq vs pr/nq and (b) v-^jriq vs KEr/riq for 
identified particle species obtained in minimum bias Au+Au 
collisions. The STAR data are from Refs. [2ll37j. 



dependent of colliding system because e is proportional 
to the pT-integrated V2 values (i.e. e = fc x W2). The 
latter proportionality has been observed for Au+Au col- 
lisions 13 nil ■ A Glauber model estimate of e [s^ gives 
/c = 3.1 ± 0.2 for the cuts employed in this analysis. This 
method of scaling leads to a scale invariant variable and 
cancels the systematic errors associated with estimates 
of the reaction plane resolution and the eccentricity. 

The resulting scaled V2 values for Cu+Cu and Au+Au 
collisions, are shown in Fig.Q^c). To facilitate later com- 
parisons with the model calculations of Ref. [l^, they 
are divided by k — 3.1. These scaled values are clearly 
independent of the colliding system size and show es- 
sentially perfect scaling for the full range of centralities 
(or e) presented. The V2 are also in accord with the 
scale invariance of perfect fluid hvdrodvnamics |23|, |27l| . 
which suggests that rapid local thermalization [9|, Il0l | is 
achieved. 

The magnitude of V2 /e depends on the sound speed 
|23l | . As a reasonable first approximation we compare our 
measured W2/e at an integrated {pt) 0.45 GeV/c and the 
results of Fig. 2 of . This results in a speed of sound 
Cs ~ 0.35 ± 0.05. Note that the calculations are done 
at fixed b=8 fm and a constant speed of sound. Thus, 
since we expect the speed of sound to vary as a function 
of time, one might view this Cs value as the approximate 
average value over the time period 2 R/cs, the time over 
which the flow develops. This value suggests an effective 
EOS, which is softer than that for the high temperature 
QGP but does not reflect a strong first order phase 
transition in which Cg = Q during an extended hadroniza- 
tion period. 

Figures El and El show that the distinctive features of 
the V2 for identified particles provide another detailed set 
of scaling tests. Fig.EJa) shows a comparison of the mea- 
sured differential anisotropy V2{pt), for several particle 



species obtained in minimum bias Au+Au collisions at 
Y^T/VTV ~ 200 GeV. The results are in good agreement 
(better than 3%) with those of our previous measure- 
ments |2J]. The values for neutral kaons (K^), lambdas 
(A) and the cascades (S) show results from the STAR 
collaboration (22^ ^7] . The STAR V2 values were multi- 
plied by the factor 1.1 to account for a small difference 
between the average centralities for minimum bias events 
from the two experiments. PHENIX and STAR V2{pt) 
results (for tt^, p(p) and K) for 10% centrality bins are 
essentially identical. 

The comparison in Fig.E^a) shows the well known par- 
ticle identification (PID) ordering of V2{pt) at both low 
and high px values. At low pT {pr ^ 2 GeV/c), one can 
see rather clear evidence for mass ordering. If this aspect 
of V2 is driven by a hydrodynamic pressure gradient, the 
prediction is that the differential V2 values observed for 
each particle species should scale with KEt- The pres- 
sure gradient that drives elliptic flow is directly linked 
to the collective kinetic energy of the emitted particles. 
For higher values of px {pt 2 — 4 GeV/c), Fig. El^a) 
indicates that mass ordering is broken and V2 is more 
strongly dependent on the quark composition of the par- 
ticles than on their mass, which has been attributed to 
the dominance of the quark coalescence mechanism for 
PT 2 - 4 GeV/c ^lillil. 

Figure [3b) shows the same V2 data presented in 
Fig. Ha) plotted as a function of KEt- Note that KEt 
is a robust scaling variable because it takes into account 
relativistic effects, which are especially important for the 
lightest particles. In contrast to the PID ordering ob- 
served in Fig.Efa), all particle species scale to a common 
set of elliptic flow values for KEt ^ 1 GeV, confirming 
the strong infiuence of hydrodynamic pressure gradients. 
For KEt ^ 1 GeV, this particle mass scaling (observed 
for all particle species) gives way to a clear splitting into 
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a meson branch (lower V2) and a baryon branch (higher 
V2)- Since both of these branches show rather good scal- 
ing separately, we interpret this as an initial hint for the 
degrees of freedom in the flowing matter at an early stage. 

Figure 121 shows the results obtained after quark num- 
ber scaling of the V2 values shown in Fig. |21 That is, V2, 
Pt and KEt are divided by the number of constituent 
quarks Uq for mesons {uq = 2) and baryons (n^ = 3). Fig. 
inia) indicates rather poor scaling for pr/nq < 1 GeV/c 
and much better scaling for pr/nq ^ 1.3 GeV/c, albeit 
with large error bars. In contrast, Fig. |3Ib) shows excel- 
lent scaling over the full range of KEx/nq values. We 
interpret this as an indication of the inherent quark-like 
degrees of freedom in the flowing matter. These degrees 
of freedom are gradually revealed as KEt increases above 
~ 1 GeV (cf. Fig.[2Ib)) and are apparently hidden by the 
strong hydrodynamic mass scaling, which predominates 
at low KEt- The fact that V2/nq shows such good scal- 
ing over the entire range of KET/riq and does not for 
PT/riq, serves to highlight the fact that hydrodynamic 
mass scaling is preserved over the domain of the linear 
increase in KEt- Fig. Otb) should serve to distinguish 
between different quark coalescence models. 

In summary, we have presented the results from 
detailed tests of hydrodynamic scaling of azimuthal 
anisotropy in Au+Au and Cu-f Cu collisions at ^Jsnn— 
200 GeV. For a broad range of centralities, eccentric- 
ity scaling is observed for charged hadrons for both the 
Cu-|-Cu and Au-t-Au systems. For a given eccentricity, V2 
is also found to be independent of colliding system size. 
The observed scaling for identified particles in Au-|-Au 
collisions, coupled with e scaling, gives strong evidence 
for hydrodynamic scaling of V2 over a broad selection 
of the elliptic flow data. For KEt ~ 1 — 4 GeV uni- 
versal hydrodynamic scaling is violated, but baryons and 
mesons are found to scale separately. Quark number scal- 
ing {v2/nq VS. KEt /riq) in this domain leads to compre- 
hensive overall scaling of the data, with substantially bet- 
ter scaling behavior than that found for V2/nq vs. pTjriq- 
The scaling with valence quark number may indicate a 
requirement of a minimum number of objects in a local- 
ized space that contain the prerequisite quantum num- 
bers of the hadron to be formed. Whether the scaling 
further indicates these degrees of freedom are present at 
the earliest time is in need of more detailed theoretical 
investigation. 
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